Low Temperature Static and Dynamic Behaviour of the Easy-Axis Heisenberg 

Antiferromagnet on the Kagome Lattice 
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The antiferromagnetic Heisenberg model with easy-axis exchange anisotropy on the Kagome lat- 
tice is studied by means of Monte Carlo simulations. From equilibrium properties, we find that 
the values of the critical exponents associated with the magnetization at the critical temperature 
Tc vary with the magnitude of the anisotropy. On the other hand, the spin-spin autocorrelation 
functions have a stretched exponential behavior with a power law divergence of the relaxation time 
at a glass-like temperature Tg ~ Tc. From non-equilibrium dynamics at a fixed temperature below 
Tg, aging effects are found which obey the same scaling laws as in spin glasses and polymers. 

PACS numbers; 75.40.Cx, 75.40.Mg 
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I. INTRODUCTION 

The Kagome spin system has attracted much inter- 
est, both theoretically and experimentally. Because the 
geometry of the lattice consists of corner sharing trian- 
gles in a layer which surround hexagons, spin systems 
are highly frustrated when antiferromagnetic interactions 
are presenti. In the case of the nearest-neighbour an- 
tiferromagnetic spin 1/2 Ising model, the ground state 
is disordered and the spin-spin correlation function de- 
cays exponentially^ at zero temperature. This behaviour 
differs from that which occurs in other periodically frus- 
trated 2d Ising lattices such as the triangular and fully- 
frustrated square lattices where the correlation function 
decays as a power Iaw2i4. Frustration also leads to a 
larger macroscopic entropy^ in the Kagome lattice com- 
pared to the triangular lattice^. 

The classical isotropic Heisenberg antiferromagnet on 
the Kagome lattice also has a ground state with macro- 
scopic degeneracy. Various perturbations such as quan- 
tum fluctuationsLSiS, or the addition of other cou- 
plings including further neighbour interactionsifi, easy 
planeii or easy-axis anisotropy^^ and Dzyaloshinski- 
Moriya interactionsi^, have a strong effect on the ground 
state manifold. There has been a great deal of contro- 
versy about whether magnetic orderiiSiiiii^xists in the 
ground state or if it remains disordered^iSiS. Two partic- 
ular coplanar and ordered states, called the q = and 
the (-s/S X -s/S) configurations, are favoured by entropy 
effects and this effect is often referred to as order by dis- 
order. The latter ordered statoiSiiLi^ is the one that 
is favoured or perhaps even a mixed disordered state of 
both structures"'-^. 

A chirality k vector characterizes these ordered con- 
figurations and is defined as the pairwise vector product 
clockwise around a triangle: 

[Si X S2 + S2 X S3 + S3 X Si] (1) 



(3\/3) 



where 1,2,3 label the three sites in the unit cell and 
form three interpenetrating triangular sublattices. The 
two ordered states mentioned above are shown in figure 
1 and correspond to configurations with uniform (ferro- 
magnetic) and staggered (antiferromagnetic) chiralities 
respectively. The structure of these ground states al- 
lows for the formation of collective zero-energy spin rear- 
rangements, called weathervane defects, that permit the 
system to explore the ground state manifold. These ex- 
citations involve the continuous rotation of the spins on 
two of the sublattices about the direction defined by the 
third. The defects may either traverse the entire lattice 
("open" spin folds) as in the q = state or form lo- 
calized loops ("closed" spin folds) as in the (\/3 x -\/3) 
stateiaJii& (see figure 1). 

It is known that most quasi-two dimensional magnetic 
materials exhibit some kind of spin anisotropy which may 
be of the easy-axis2£ or easy-plane^^ type. Easy-plane 2d 
magnets have attracted attention due to the possibility of 
a topological Kosterlitz-Thouless phase transition which 
may exhibit glassy behavior different from that found in 
conventional site-disordered systems^. The amount of 
interest devoted to easy-axis magnetic systems has been 
considerably smaller, especially with regard to the study 
of dynamical properties. Kuroda and Miyashita™ (KM) 
have previously studied the Ising-like Heisenberg anti- 
ferromagnet on the Kagome lattice using Monte Carlo 
methods. They have shown the existence of a phase tran- 
sition at very low temperature with an exotic ordered 
phase which has no spatial long-ranged order and hence 
shares some similarities to spin glasses. In the present 
work we study both the equilibrium and dynamic prop- 
erties of this model. 

The paper is organized as follows: in section II we 
describe the model and our methods of calculation. Sec- 
tion III describes the equilibrium properties of the model 
and the critical exponents deduced from finite size scal- 
ing. Results for the spin-spin autocorrelation function 
are also given in this section. A good fit is obtained us- 
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(a) 



(b) 




FIG. 1: Isotropic Heisenberg antiferromagnet on the Kagome 
lattice: (a) the q=0 state in which the spins on each of the 
sublattices are parallel to each other and make an angle of 
120° with the spins on the other two sublattices and (b) the 
(•s/3 X y/S) structure has a larger unit cell. The + and — on 
the triangles indicate the chirality and the circles describe the 
open and closed spin folds in the states q=0 and (\/3 x \/3) 
respectively. 

ing a functional form which has been used to describe 
3d spin glasses^^ and other complex systems^Si^iS^. A 
characteristic timescale r diverges with a power law at 
a temperature close to the static critical temperature. 
In section IV, we study the dynamical behavior below 
the critical temperature. The spin-spin autocorrelation 
function exhibits aging effects characteristic of glasses. 
We summarize our results in section V. 



II. MODEL AND METHODS 

The model is described by the following Hamiltonian 

H^jY,{s^s^ + sfsy + As^s;). (2) 

where {Sf,a = x,y,z) represents a classical three com- 
ponent spin of unit magnitude located at each site i of 
a Kagome lattice and the exchange interactions are re- 
stricted to nearest-neighbour pairs of sites. The parame- 
ter A describes the strength of the exchange anisotropy. 
We restrict our attention to the case where A > 1 rep- 
resents an easy-axis anisotropy. The limit A — s- 1 cor- 



responds to the isotropic Heisenberg model whereas the 
limit A —^ oo corresponds to an infinite spin Ising model. 
The model has a macroscopic ground state degeneracy 
for all ^ > 1 with a ground state energy per site given 

3 A+1 ■ 

The ground state of the system for A > 1 corresponds 
to a configuration in which the spins on each triangle 
form a distorted 120° planar state with a net nonzero 
magnetization in the zb^-direction whose magnitude is 
related to A as |m^| = (^^). The local chirality n is nor- 
mal to the plane of each triangle but it does not show any 
evidence of long range order. Miyashita and Kawamura— 
(MK) have previously studied the same model on the tri- 
angular lattice and observed a non-trivial degeneracy re- 
lated to the rotation of the magnetization vector in the 
plane of the triangles. In contrast to the corner-shared 
triangles of the Kagome lattice, the triangular lattice 
shares edges and there is a q = sublattice order at 
T = 0. There are also two sequential finite temperature 
topological phase transitionsSLSS with the lower transi- 
tion corresponding to the onset of a power law decay of 
the chirality. The same degeneracy arguments of MK 
apply to the Kagome case but at finite temperature this 
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FIG. 2: A Monte Carlo snapshot of the spin configurations 
at low T for the case A = 2. The lower panel shows the 
X — y spin plane. On each triangle there is one spin in the 
z-direction but there is no spatial sublattice order. The upper 
panel shows the total magnetization on each upward triangle 
for the same spin configuration. There is a net magnetization 
in the z-direction. 
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degeneracy seems to be lifted by an order from disorder 
effect and the z-axis is preferred. This could be related to 
the fact that, in the Heisenberg limit A ^ 1, the weath- 
ervane modes described in figure 1 can rotate about any 
axis whereas for A > 1 these excitations select the z- 
axis. Monte Carlo snapshots of the spin configurations 
at very low temperatures reveal an exotic phase for which 
there is no evidence of long-ranged spatial order of the 
individual spinsi^. Both weathervane modes and other 
localized defects are observed in these snapshots. Figure 
2 shows a typical configuration at very low T in the case 
A = 2. The lower panel shows the x — y components of 
the individual spins and the upper panel shows the cor- 
responding z component of the magnetization on each 
upward traingle. A localized defect is observable which 
corresponds to a triangle with two spins up and one down 
and zero chirality. Both the magnetization and heat ca- 
pacity exhibit critical behaviour at a finite temperature 
Tj, corresponding to the breaking of the z-axis up-down 
symmetry and is similar to the 2d ferromagnetic Ising 
model. 

We employ Monte Carlo methods using a single spin 
flip heat bath algorithm to study lattices containing N 
spins with periodic boundary conditions. The number of 
spins is related to the number of unit cells as = 3L^, 
where L is the number of up-triangles in the horizontal 
direction. We have calculated various thermodynamic 
observables such as the internal energy, the specific heat, 
the z-component of the magnetization as well as the asso- 
ciated susceptibility and Binder cumulant^^. Our numer- 
ical data are analyzed by using finite size scaling theory 
and the histogram methodii^Si^ to extract the critical 
exponents of this model. A reweighting method is com- 
bined with our single spin flip algorithm in order to ob- 
tain the observables as continuous functions of tempera- 
ture near Tc- We first measure the specific heat using the 
MC simulations on a discrete temperature grid and this 
step yields an estimate of the temperature Tq at which 
the specific heat is maximum. Using this estimate, the 
histogram ^^{E) of the number of spin states with en- 
ergy E is constructed from MC runs at the temperature 
To over a large time interval At. This procedure allows 
us to obtain the average value of any observable Q as a 
continuous function of temperature T near Tq as follows 



(3) 



where Q{E) is the microcanonical average of the ob- 
servable. This method has been used quite succes- 
fuUy to extract critical exponents of both discrete^&and 
continuousiiSiiS^ spin models. 

We have also studied dynamical properties of this 
model by considering the double-time spin-spin autocor- 
relation function 
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FIG. 3: (a) The order-parameter Binder cumulant Ul (12 < 
L < 60) plotted vs T obtained by optimized reweighting in the 
case A=2. (b) Estimates for Tc plotted vs inverse logarithm 
of the scale factor b — L' /L. 



To measure this quantity at a given temperature T we 
start from a random configuration at high temperature 
and rapidly quench to the working temperature T. We 
then wait for a time and measure the autocorrelation 
function C(t, tw) for subsequent times t. The results are 
averaged over many random initial states. In equilibrium 
one expects C{t, t^) to be independent of and it is only 
in equilibrium that one can define a meaningful timescale 
associated with relaxation. In the aging regime, C{t,tw) 
is waiting time dependent^S*^. Aging is a general phe- 
nomenon which occurs in a wide variety of off-equilibrium 
materials, as for example glasses. The phenomenon has 
been widely studied in disordered systems such as spin 
glasses^^., frustrated systems^ and in the phase ordering 
kinetics of the Ising ferromagnetS^, and is associated with 
a slow domain dynamics. 



III. EQUILIBRIUM PROPERTIES 



Cit,t^) = ^ < E + > (4) 



We have used our MC method to study lattice sizes 
L = 6,12,18,24,36,48,60 and have used 1 - 5 x 10^ 
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FIG. 4: Finite size scaling dependence of the critical properties for A = 2 (a) The order parameter Mz at Tc, (b) the temperature 
derivative of the Binder cumulant associated with M^, at Tc, (c) the maximum of the susceptibility y^"-^ and (d) the specific 
heat maximum (7™-"^ plotted as a function of L on a log-log scale. 



Monte Carlo steps (MCS) for performing our measure- 
ments after discarding the first 5 x 10^ MCS to reach ther- 
mal equilibrium. In the reweighting analysis, it is impor- 
tant to take At as large as possible to have good statistics. 
We have used Ai = 1.2 x 10*^ to At = 2.6 x 10^ MCS for 
small and large lattice sizes respectively. Also, because 
the energy is continuous, we have used both 10000 and 
30000 bins for the histograms in order to check that the 
size of the bins did not affect our numerical results. In 
the results that follow, the magnitude of the exchange 
constant J is set equal to unity. 

The critical temperature, Tc, can be determined by 
comparing the reduced Binder cumulant of the magne- 
tization, Ul = 1 - < >/3< Ml >2, for lattices of 
size L with lattices of size L' = hh as shown in figure 
3(a) for the case of the anisotropy parameter A = 2. 
In the limit of large system sizes, the cumulants should 
cross at the critical temperatur o^^i'^^ and have a com- 
mon value Ul = U* . However, due to finite size effects, 
it is necessary to extrapolate the crossing points to the 
limit b — > 00^. Our results using the Binder cumulant 
crossing methodSS to estimate the critical temperature 
are presented in figure 3(b). The points represent the 
temperatures at which the order parameter cumulant for 



L' crosses the cumulant for T = 12 or T = 15. There is 
considerable scatter in the data and care must be taken 
to use only results with L' sufficiently large to be in the 
asymptotic region where a linear extrapolation is justified 
(i^ < 2.2). Using this method, the critical temperature 
is"estimated to be T^ = 0.077 ± 0.001. This value is 
slightly lower than the one obtained by KM from phe- 
nomenological renormalization of the magnetization but 
lies within their error bars. No hysteresis is observed in 
the order parameter nor in the energy near the critical 
region. In addition, no double peak structure was found 
in the energy histograms and the Binder energy cumu- 
lant evaluated at Tc yielded the result U* = 0.666665(7) 
for large L, consistent with the value | expected for a 
continuous transition"^^. 

Finite-size scaling results for the order parameter M^, 
the first temperature derivative of its Binder cumulant, 
the susceptibility x = ^(< > - < and 
the specific heat are shown in figures 4 (a-d) respec- 
tively on a log-log scale. According to the standard the- 
ory of finite size scaling, the equilibrium magnetization 
Mz should obey the relation M^ ^ L^^^/" for sufficiently 
large L. Figure 4(a) shows our results of a finite size scal- 
ing analysis for the order parameter M^. Excluding the 
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TABLE I: Results for the static critical temperature Tc and 
the exponents 'y/u, fijv and v for various values of the 
anisotropy A. 



A 


Tc 


7/^ 


(ijv 


Ijv 


1.1 


0.036 ± 0.006 


1.44 ± 0.07 


0.30 ± 0.05 


0.97 ±0.06 


1.5 


0.067 ±0.001 


1.61 ±0.05 


0.22 ± 0.04 


1.03 ±0.04 


2 


0.077 ±0.001 


1.64 ±0.03 


0.20 ± 0.02 


0.98 ±0.02 


3 


0.076 ±0.001 


1.66 ±0.05 


0.18 ±0.03 


0.99 ±0.03 


5 


0.064 ± 0.002 


1.67 ±0.04 


0.17 ±0.03 


1.01 ±0.03 


8 


0.052 ± 0.003 


1.67 ±0.04 


0.16 ±0.01 


1.09 ±0.06 


14 


0.037 ±0.001 


1.70 ±0.05 


0.14 ±0.01 


1.02 ± 0.04 


20 


0.030 ±0.001 


1.73 ±0.03 


0.13 ±0.01 


1.03 ± 0.04 


30 


0.022 ± 0.002 


1.72 ±0.04 


0.13 ±0.02 


1.02 ± 0.04 




smallest two lattice sizes L = 12 and 15 from the fitting 
procedure, we obtained the value of the exponent ratio 
P/v = 0.20 ± 0.02 which is significantly larger than the 
2d Ising value P/v = 1/8 . 

The behavior of the reduced Binder cumulant Ul at 
the critical point can be used to find the value of the 
critical exponent v. Finite size scaling theory predicts 
at Tc that Ul = UoitL^/") with t =\ 1 - T/T^ \ and 
the temperature derivative of Ul at should obey the 
relation U'^{Tc) = L^/^U^iQ). In figure 4(b) we show 
that this prediction is obeyed quite well. The value of 
the static exponent v obtained using a least-squares fit is 
1/v = 0.98 ± 0.02, which is remarkably close to the two 
dimensional Ising value v = 1. The magnetic suscepti- 
bility X has the scaling form x ^ L^^" and figure 4(c) 
shows a least squares fit to our results using this form 
and we find 7/1/ ~ 1.64 ± 0.03 which is smaller than the 
2d Ising value 7/1^ — 7/4 and also the value obtained by 
KM. The specific heat was also calculated but it is much 
more difficult to analyze because of the small number of 
points used and the scatter in the data was too large to 
extract a reliable estimate for a/v . Cmax should scale in 
the critical region as CJ^""^ ~ Co + aL""^" , with Co rep- 
resenting the regular part. From our fitting procedure, 
C'o = yields the best straight line for large sizes with 
slope a/v = 0.07 ±0.04 . 

The same analysis has been carried out for other val- 
ues of the anisotropy A and the results are summarized 
in Table I. The critical temperature increases from zero 
for small deviations of A from unity, attaining a maxi- 
mum near A 2 and decreasing to very low temperatures 
for large A in agreement with KM. In the limit A ^ 00, 
the transition temperature approaches zero as in the case 
of the spin 1/2 antiferromagnetic Ising model. It is not 
clear whether or not there is an ordered ground state 
in the limit A ^ 00. The critical exponents are plot- 
ted as a function of the anisotropy parameter A in figure 
5. The values of the critical exponents P/v and ^/v ap- 
pear to depend on the value of A and approach the usual 
2d ferromagnetic Ising values as A becomes very large. 



FIG. 5: Variation of the critical exponents with the anisotropy 
A obtained from the equilibrium properties. The symbols 
represent the measured values and the lines indicate the values 
expected for the 2d Ising model. 

Both universality and weak universality^^ are violated in 
this system. This non-universal behaviour of Ising-like 
exponents has also been reported for a two dimensional 
system of XY spins interacting via both ferromagnetic 
and antiferromagnetic bonds in the presence of an ap- 
plied magnetic field which reduces the symmetry 0(2) in 
spin space to Z-^. 

In order to understand the nature of this exotic struc- 
ture in more detail, we have also looked at the spin-spin 
autocorrelation function using equilibrium dynamics at 
high temperatures where C(i, i^) becomes independent 
of for > 10"*. As the temperature is reduced, the 
relaxation of the spins becomes slower and deviates from 
a simple exponential form. We fit our data with the fol- 
lowing functions^, 

/oW - ^^^p(-(Vt)''^^^) (5) 

where a,b{T),T{T),x{T) are all fitting parameters. As 
shown in figure 6(a) for A — 2, our data are fairly well 
described by this functional form and we are able to ex- 
tract the temperature dependent relaxation time t{T) 
and the exponent b{T). The behavior of these param- 
eters is depicted in figure 6(b). The relaxation time r 
appears to increase sharply as T is reduced and can be 
fit quite accurately by a power law divergence of the form 

T~(r-T,)-'^ (6) 

The glass temperature, Tg, identified from the power law 
divergence of the relaxation time is slightly lower than 
the critical temperature Tc obtained from the equilibrium 
measurements for values of A which are not too large. For 
A ~ 2, it seems that there is a decoupling between the 
local spin degrees of freedom and the net magnetization 



6 



C(t) 0.1 : 



0.01 

9000 r- 

8000 - 

7000 - 

6000 - 

5000 - 

4000 - 

3000 - 

2000 - 

1000 - 

- 




100 1000 
t (mcs) 



10000 100000 



(b) 



0.15 0.3 0.35 0.3 0.35 



0.05 



0.1 0.15 0.2 0.25 0.3 0.35 0.4 

T 



FIG. 6: (a) The spin-spin autocorrelation function for various 
temperatures, T = 0.4(+), T = 0.15(*), T = 0.09(x),r = 
0.075(o) for A = 2. The lines are fits to form given in (5) (b) 
The associated relaxation time as a function the temperature 
together with a fit to the power-law form (6). The inset shows 
the temperature dependence of the exponent 6(T) . 



on each triangle. The composite spin variable orders at 
Tc and the local spin variables enter a glassy phase at 
Tg < Tc- At larger values of A the glass temperature Tg 
approaches the critical temperature Tc- As will be shown 
in the next section, Tg signifies the onset of aging phe- 
nomena and it is impossible to define a timescale below 
this temperature since the system is in a frozen glassy 
state. 

The exponent 6(T) is plotted in the inset of figure 6(b). 
It is temperature dependent and lies in the range bi < 
b < 1 ioT Tg < T < T* where 61 = 0.75, = 0.071 
and T* = 0.1 for ^ = 2 . The non-exponential behavior 
sets in at temperatures below T* . The parameter x{T) 
which characterizes the short time behaviour lies in the 
range < x{T) < 0.1 and decreases with temperature. 
A summary of our relaxation results for a few values of 
A are given in Table II. 

The relaxation time exponent is also nonuniversal 
and decreases in value for larger values of A. Both non- 
exponential relaxation and a diverging relaxation time 



TABLE II: Results obtained from high temperature equilib- 
rium dynamics of the spin-spin correlation function C{t,tw): 
Tc is the transition temperature obtained from the equilib- 
rium properties, T* is the temperature at which C{t,tm) first 
has a non-exponential behaviour; Tg is the critical tempera- 
ture obtained from (6) where the relaxation time diverges; 61 
is the lowest value of the exponent b{T) in (5) near Tg and (j> 
is the relaxation time exponent in (6). 



A 




To 


T* 


bi 





1.1 


0.036 ±0.006 


0.032 ± 0.001 


0.050 


0.76 


1.5 


2 


0.077 ±0.001 


0.071 ± 0.002 


0.100 


0.75 


1.3 


8 


0.052 ±0.003 


0.05 ±0.001 


0.055 


0.85 


0.74 


30 


0.022 ±0.002 


0.022 ± 0.001 


0.026 


0.86 


0.55 



are features of glasses^Si and this behaviour has previously 
been seen in frustrated systems without disorder^'^. 



IV. OFF-EQUILIBRIUM DYNAMICS 

In order to further study the dynamics of this model, 
we have carried out some numerical experiments focused 
on revealing the presence of slow dynamics in conjunc- 
tion with history-dependent phenomena which is gener- 
ally referred to as 'aging'. These features are most eas- 
ily found in simulations of the two-time autocorrelation 
function C(t,tu,) which shows an explicit dependence on 
both times t, over a wide range of time scales. Aging 
can be observed in real systems through different exper- 
iments. A typical example is the zero-field cooling ex- 
periments in which the sample is cooled in zero field to 
a subcritical temperature at time t = Q. After a waiting 
time tw a small magnetic field is applied and subsquently 
the time evolution of the magnetization is recorded. It is 
often observed that the relaxation becomes slower as the 
waiting time tw is increased. 

We have measured the behavior of C{t, i^) as a func- 
tion of the observation time t, for different values oit^^A 
and T. We have used 1.5 x 10^ MCS with a lattice size 
L = 36 and averaged the results over approximatelylOG 
different trials. At high temperatures T > Tg, we found 
that the system does not exhibit aging since for any value 
of A the autocorrelation functionC(t, t^) is homogeneous 
in time and independent of t^, . 

In figures 7(a-b) , the behaviour of the autocorrela- 
tion function clearly confirms the presence of aging in 
this model for all values of A > 1 at very low temper- 
atures T < Tg. For large waiting times and t << t^, 
the correlations are independent of the waiting time t^. 
However, for t > tw, the curves show an explicit depen- 
dence on both times indicating that equilibrium has not 
been attained within the time of the simulation and the 
correlation falls to for t ^ 00. This scenario has been 
called weak ergodicity breaking^Mi. The fiuctuation dis- 
sipation theorem holds for short times but is violated at 
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longer times. 

We have attempted to find an appropriate scaling law 
for the aging curves. A knowledge of the scaling form 
could give some insight into the nature of the underlying 
dynamical process, even if there is no theoretical basis for 
determining the scaling functions. The simplest scenario 
is naive aging of the form 

C(t,i.) = /(f) (7) 

in the region where both t and are large. 

In figure 8 we show C{t,tw) as a function of t/t^ to 
see if this naive form of scaling holds. Except for the 
largest value of = 50000, we observe a departure 




0.2 ^ — — — — "J— — "J— — — — ^ 
1e-05 0.0001 0.001 0.01 0.1 1 10 100 1000 

K 

FIG. 8: Autocorrelation function C{t, tw) vs t/t^ at T = 0.04, 
for A = 2. 



from naive scaling with the function C{t,t^) increasing 
when t > tw and decreasing for t < a,t fixed val- 
ues oit/tyj as increases. A 'superaging' behaviour as 
observed in mean field spin glasses or the Sherrington 
Kirkpatrick (SK) model''^ could be expected for smaller 
waiting times but the fact that the curve for the largest 
waiting time, t^] = 50000, lies below the next smaller 
value, tw = 25000, could be explained as follows: the 
system is in a 'subaging' scaling region where the relax- 
ation of older systems becomes faster when plotted versus 
t/tw although, when plotted versus t, the older the sys- 
tem appears to exhibit slower relaxation. Indeed, as seen 
in figure 9(a), a good collapse in the asymptotic region of 
the largest waiting times is obtained by using a variable 
used in glassy polymers^ and recently in a topological 
spin glass^i*^, namely [{t -\- twY~^ — — /i) where 

in our case the value /i = 0.8 < 1 is used. This subaging 
effect has not only been observed in glassy polymers^^ 
but also in 2d site-disordered spin glassed. 

For larger values of the anisotropy, our analysis of the 
largest waiting times shows that simple scaling holds, 
compatible with the full aging scenario. Figure 9(b) 
shows the results for A — 30. This simple scaling has 
also been observed in the 3d Edwards- Anderson spin 
glasa4",4ii45. 



V. SUMMARY 

In this work we have performed a numerical study of 
the two dimensional easy-axis Heisenberg antiferromag- 
nct on the Kagome lattice by computing its static and 
dynamic properties. From the static properties, we have 
extracted the critical temperature and the critical expo- 
nents associated with the magnetization, the susceptibil- 
ity and the correlation length respectively. Our result for 
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function, C{t,tw), at high and low temperatures respec- 
tively. We have found that the high temperature equi- 
librium correlation function is described very well by the 
function at~^exp{—{t/T)^) suggested by Ogielski2i over 
the entire time and temperature range. Non-exponential 
relaxation sets in at a temperature T* > Tc- The re- 
laxation time, T, increases according to a power law and 
diverges at a temperature Tg < Tc where a transition to 
a glassy phase is located. Figure 10 shows our results ob- 
tained for the temperatures Tc, Tg and T* obtained from 
the statics and dynamics plotted as a function of A. KM 
also identified a broad maximum in the specific heat for 
A > 2 at a lower temperature than those in figure 10 
which they attribute to a local degree of freedeom since 
the peak does not exhibit any size dependence. This lo- 
cal degree of freedom could be the weathervane mode or 
the defect observed in figure 2. 

Maegawa et. al.^^ have reported an observation 
of successive phase transitions in the Kagome systems 
RFez{OH)Q[SOi)2[R = NHi,Na,K]. Susceptibility 
cusps are observed at two closely spaced temperatures 
which are about 10% of the corresponding Curie Weiss 
temperatures. These transitions could be explained in 
terms of the ordering of the magnetization on each tri- 
angle at the upper temperature Tc followed by a spin 
freezing of the local spins at the lower temperature Tg. 

Below Tg, we have found clear evidence for the pres- 
ence of aging effects in the autocorrelation function from 
off-equilibrium dynamics. The spin-spin autocorrelation 
function depends on both times and the dynamics be- 
comes slower for larger waiting times. An analysis of 
the autocorrelation functions from scaling forms used in 
polymer glasses and spin glasses has shown different be- 
haviour. Namely, a sub-aging behavior at low values of 
A is seen where the relaxation time of the system grows 



FIG. 9: Data collapse of the curves shown in figures 7(a,b) at 

large waiting times for (a) T = 0.04 and yl = 2 as a function 0.12 i • ' ' — • — • — ■ ■ . • — — 

of the time reduced variable [{t + t^uY''^ — tir'']/(l — ^) with j 
/.i = 0.8 and (b) T = 0.01 and A = 30 as a function of t/t™. Tg 

0.1 - 

the critical temperature obtained from the Binder cumu- 08 - j, 

lant method is in agreement with that obtained by KMi^ * 
within statistical errors. On the other hand, our results 
for the critical exponents indicate that this system has 
a nonuniversal phase transition. Namely, the values of 
the exponents (3 and 7 associated with the order param- 
eter and the susceptibility vary with the magnitude of 
the easy-axis anisotropy A. However, a and v remain 
unchanged and correspond within errors to the 2d Ising 
values and thus weak universality is also violated. Al- 
though the magnetization indicates a finite Tc, Monte 
Carlo snapshots of the individual spins below this tem- ^ 
perature do not indicate any long ranged spatial order. 

Rather, the individual spins appear to be in a frozen state FIG. 10: The temperatures T*, and Tg plotted as a func- 
similar to a glass. tion of the anisotropy A using the values listed in Tables I 

We have studied the two-time spin-spin correlation ^^'^ lb 
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more slowly than the waiting time tw as observed in 2d 
spin glasses^, polymer glasses^ and structural glasa^, 
whereas for large values of A, a full aging behavior de- 
scribes the data well, where the relaxation time of the 
system scales with its age as observed in 3d spin 
glasses^S^. We are currently extending our investiga- 
tions of this system to include the effects of a small ap- 
plied magnetic field on the aging behaviour. This will 
enable us to check whether the fluctuation dissipation 
theorem is violated and to study the long-term memory 
of this model. 
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